ORIGINS OF DIFFUSION 

ANTTI KUPIAINENi 

Abstract. We consider a dynamical system consisting of subsystems indexed by a 
lattice. Each subsystem has one conserved degree of freedom ("energy") the rest being 
uniformly hyperbolic. The subsystems are weakly coupled together so that the sum 
of the subsystem energies remains conserved. We prove that the long time dynamics 
of the subsystem energies is diffusive. 

o": 

o 

CN ; 1. Diffusion from conservative dynamics 

>>■ 

c^ ' One of the fundamental problems in deterministic dynamics is to understand the 

microscopic origin of dissipation and diffusion. On a microscopic level a physical sys- 
tem such as a fluid or a crystal can be described by a Schrodinger or a Hamiltonian 
dynamical system with a macroscopic number of degrees of freedom. Although the 
' c-i '' microscopic dynamics is reversible in time one expects dissipation to emerge in large 

P^. spatial and temporal scales e.g. in the form of diffusion of heat or concentration of 

r-| ! particles. 

^y ' To fix ideas, consider a Hamiltonian dynamical system i.e. a Hamiltonian flow on a 

symplectic manifold M. For the present purpose it suffices to consider M = M^" with 
position and momentum coordianates q,p ^ M". The Hamiltonian flow (pt G DiffM 
generated by the vector field {dpH, —dqH) where if : M — )■ M is the Hamiltonian or 

>• ' energy function preserves the energy 

oo ; 

OO . Ho(l)t = H 

lO ■ 

(^ ■ i.e. the flow preserves the constant energy sets Me = {{(l,p) '■ H{q,p) = E}. 

if^ . On the other hand, the simplest diffusion process is given by the heat equation 

§1 (1) dtE{t,x) = KAE{t,x) 

and the associated semigroup -ipt = e*^*^. Unlike for the reversible (pt where (p-t = 0^"^, 
ipt has no inverse and describes dissipation. Physically, the energy function E{t, x) 
^ ■ describes a macroscopic energy density i.e. a coarse grained function of microscopic 

dynamical variables, the positions and momenta of the underlying Hamiltonian dy- 
namics. The question we wish to pose is how does this dissipative dynamics tpt arise 
from the conservative one (pf 

A concrete physical system where diffusion occurs is a fluid. In classical mechanics 
this is microscopically modeled by a Hamiltonian system whose flow gives the trajec- 
tories of the fluid particles {qi{t) , Pi{t)) G M^ x M^, i = 1 . . . A^. A typical Hamiltonian 
function is given by 

(2) Hii.P) = T.§^ + T.^(^^-^^) 

i ij 



^ 
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consisting of the kinetic energy of the particles of mass m and a pair potential energy 
of interaction of the particles. Let the energy of the i:th particle be defined as 

(3) ^» = || + IE^(^^-^.) 

so that H = X]j^j- W^ can then define the energy density as the distribution 

(4) E{t,x) = J2('iK{t){x) 

i 

where Sg is the Dirac mass at q. Since J E{t,x)dx = Yli^i — H and H = one 
concludes 

(5) E{t,x) = V ■Jit,x) 

for a certain distribution, the energy current, depending on q{t),p{t). Eq. ([5]) is a local 
conservation law deduced from the global energy conservation. In the case of the fluid, 
there are two other similar local conservation laws related to global momentum and 
particle number conservation laws. This leads to a richer set of macroscopic laws in the 
case of the fluid than the diffusion equation for the energy (in particular these include 
the Navier-Stokes equations). 

2. Coupled oscillators 

Thus, to understand the origins of diffusion one should look for systems with just 
one local conservation law eq. ([5]). There has been a lot of work in recent years around 
these questions in the context of coupled dynamics i.e. dynamical systems consisting 
of elementary systems indexed by a rf-dimensional lattice Z'^. The total energy E of 
the system is a sum ^^ E;j. of energies E^ which involve the dynamical variables of the 
system at lattice site x and nearby sites. The physical situation to keep in mind is then 
thermal conduction in a crystal lattice (i.e. a solid). 

Two types of systems have been considered. In the first case at each lattice site 
we have an oscillator and the oscillators at neighboring sites are coupled together. 
Typically one considers the system where the forces are weakly anharmonic. In the 
second case at each lattice site one puts a chaotic system and weakly couples the 
neighboring systems. Let us start with the former case. 

The setup resembles that of the fluid above, but now the "particle" positions q^ are 
indexed by the lattice, x G A C Z*^ where A is a finite subset, say a cube, and they 
describe the deviation of an atom from its equilibrium position at x. A simple classical 
mechanical model for this is a system of coupled oscillators 

(6) ^A(g,p) = E(^ + ^(^-))+ E ^fe-^.) 

xGA \x-y\=l 

where U is a pinning potential which we assume tending to infinity as |g| — )■ oo. The 
potential V describes interaction of the atoms in nearest neighbor lattice sites and is 
taken attractive. A challenging model is obtained already by taking 

(7) Viq) = q', U{q) = q' + X\q\' 
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and further simplifying by taking g^, G M instead of W^. Then a lattice version of eq. 
(|5]) holds with the current given by 

(8) J,,{x) = -^{p^+f,+p^)V'{q^+^ - q^). 

In what sense should we expect the conservative dynamics ([5]) give rise to a diffusive 
one as in eq. ([T])? The answer is that this should happen for typical initial conditions 
(g(0),p(0)) G Ma with respect to a specific measure on the phase space M, 



A 



b2|A| 



and under a proper scaling limit. 

Recall first that the Hamiltonian dynamics preserves the Lebesgue measure m\ on 
M\. Since also H\ is preserved so is the Gibbs measure (or equilibrium measure) 

where /9 > as well as its (thermodynamic) limit /i/j = livuf^^id iij^t^. Let us now replace 
the (inverse) temperature parameter /3 by a spatially varying one. Let h G C^(M'^) and 
P > 1 1^1 loo- Write as in the fluid case 



ffA = E' 



xeA 
, Cx describing the energy contributed by the oscillator at x. Pick a scaling parameter 
L eN and set (3l{x) = (3 + b(x/L). Let fi^^' be the thermodynamic limit of the measure 

Construction of this limit poses no problems if A > in eq. (jTj) is small enough, fx^^^ 
is not invariant under the dynamics which maps it to fil = fj.^^'^ o 0^^. However, one 
expects that as t — )■ oo there is return to equilibrium i.e. /ij — > /i^. The diffusion 
equation is expected to govern this process in the following sense. 
Let / G C^(M+ X W^) and consider the random variables 

(9) ei(/) = L-'^-2 J2 f{t/L-\x/L)eMt).P{t))- 

(t,x)gZ+xZ<< 

The statement of the hydrodynamic limit is then: with probability one in the sequence 
of measures fj.^^\ ^lif) converges to J fit, x)E{t, x)dtdx where E is the solution to the 
nonlinear diffusion equation 

(10) dtE = V ■ {k{E)VE) 

where n{E) is a smooth positive function. The initial condition E{0, ■) is determined 
by the function b. Thus upon coarse graining and scaling the equation ([5]) turns to eq. 
flTOj) . almost surely in the initial conditions of the underlying microscopic variables. 

The proof of the hydrodynamic limit in our model is beyond present mathematical 
techniques. The existing techniques require the presence of plenty of noise in the 
system. A simpler problem would be to establish the kinetic limit. This is a weak 
anharmonicity limit. We replace A in eq. ([7]) by X/yL and and consider the measures 
^\^^ . As L — )■ oo we expect these measures to become gaussian whose covariance upon 
spatial scaling satisfies a Boltzman equation. More precisely, denote {qx,Px) by (p{x). 
Then it is conjectured that 

(11) lim / 0(Lx + y)(piLx - y)fif^\d(f)) = G{t, x, y) 
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exists and the Fourier transform of G{t,x,y) in y, G{t,x,k) satisfies the so called 
phonon Boltzman equation 

(12) dtG{t, X, k) + Vu{k) ■ VG{t, x, k) = I{G{t, x, ■)) 

where / is a nonlinear integral operator and u:{kY is the Fourier transform of the lattice 
operator 2(— A + 1), see [T]. Proof of these statements is still open and a considerable 
challenge (for some progress see [2]). Derivation of a hydrodynamic equation of the 
type fITOl) from the Boltzman equation flT2l) has been carried out [3] , see also [1] where 
an attempt to go beyond the kinetic limit was carried out. 

3. Coupled chaotic flows 

A second class of models deals with a complementary situation of weakly coupled 
chaotic systems [5], [6], [7]. The setup is as follows. Let {M,H) be a Hamiltonian 
system i.e. M is a symplectic manifold and if : M — t- R. Let, for each x G Z'^ (M^, H^) 
be a copy of {M,H). Let /i : M x M — )■ R and for each x,y e 11^, \x — y\ = 1 let 
h^y : Mj X Mj^ -)■ R be a copy of h. Let A C Z'^ be finite and Ma = x^^z^M^. The 
coupled flow is the one on M\ generated by the Hamiltonian 

(13) Ha = J2h.+ Y1 ^^-y 

X'gA \x—y\ = l 

Of course, the coupled oscillators of the previous section are of this form. There, 
the system (M, H) is integrable, and the diffusive dynamics is the consequence of 
coupling and anharmonicity. In the present discussion we wish to take (M, H) chaotic. 
Examples are Anosov systems or billiard systems. E.g. in the former case the flow on 
M generated by H has diniM — 2 non-zero Lyapunov exponents and two vanishing 
ones corresponding to the Hamiltonian vector field and Vif . 

When the coupling parameter A is zero {M\,Hj>^) has 2|A| vanishing Lyapunov ex- 
ponents. For A 7^ one expects that for a large class of perturbations h the only 
constant of motion is i^A and the system has only two vanishing exponents. However, 
zero should be near degenerate for the Lyapunov spectrum and these long time scale 
motions should be at the origin to diffusion in the A — t- Z"^ limit. 

Rigorous results on such Hamiltonian systems are rare: in ^ ergodicity is proved in 
a one dimensional model. However, it seems very difficult to get hold of the Lyapunov 
spectrum and it is far from obvious how such knowledge would turn into a proof of 
diffusion in these systems. I want to argue that a more fruitful approach is to study the 
local energy conservation law ([5]) and try to show that the chaotic degerees of freedom 
act there like a noise that redistributes locally the energy. To probe such an idea it is 
useful to turn to a discrete time version of our model i.e. to study iteration of a map 
rather than a flow. 

4. Coupled chaotic maps 

A discrete time version of the coupled flow setup of the previous section is called a 
Coupled Map Lattice (CML). Now the local dynamical system is a pair (M, /) where 
M is a manifold and / : M — )■ M. Again for each x G Z'' (M^, /^) is a copy of (M, /) 
and {Mj>^, f/<^) with f\ = XxgaIx is the product dynamics. The CML dynamics is a 
suitable local perturbation of the product dynamics. 
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Our choice of M and / is motivated by the coupled chaotic flows discussed before. 
A discrete time version (say given by a Poincare map) of a billiard or Anosov flow 
has one vanishing Lyapunov exponent corresponding to the conserved energy and the 
remaining ones nonzero. We model such a situation by taking for the local dynamics 
the manifold of form M = ]R_|_ x A^ with A^ another manifold. Let us denote the 
variables at the lattice site x E Z'^ hj {E{x), 0{x)) e IR+ x A^. We call the non- negative 
variables E energy and postulate them to be conserved under the local dynamics: 

(14) {E{x), eix)) -> iEix),gi9ix), E(x))) 

for each x G Z"'. 

6 E N are the fast, chaotic variables. In the billiard case the dynamical system 6 — t- 
g{6, E) is uniformly hyperbolic for any fixed E. We will model this situation by taking 
g{6, E) = g{6) a fixed chaotic map, independent of E. Examples are A^ = T^ = M/Z 
and g an expansive circle map, e.g g{d) = 26 and A^ = T^ = M^/Z^ and g a hyperbolic 
toral automorphism. 

We should stress that the E independence is the most serious simplification in this 
setup. In a realistic Hamiltonian system, such as the billiards the E dependence of g 
can not be ignored. Indeed, it is obvious that as i? — )■ the Lyapunov exponents of 
g{-,E) also tend to zero since E sets the time scale. 

The CML dynamics is a perturbation of the local dynamics (fT^ . Let us use the 
same notation {E, 6*) G Ma = M^ x A^^. Then F : Ma -^ Ma is written as 

(15) F{x, E, 6) = {E{x) + fix, E, 9),g{9{x)) + h{x, 9)). 

Here / and h are small local functions of {E, 9) i.e. they depend weakly on {E{y), 9{y)) 
for |x — yl large as we will specify later. 

/ is however constrained by the requirement that the total energy ^^E{x) is con- 
served. This follows if 

5^/(x,E,^)=0 

X 

for all E,9. A natural way to guarantee this is to consider a "vector field" J(x) = 
{J^(x)}^=i^...^d and take 

(16) /(x, E, 9) = {V- J)(x, E, 9) := 5^( J'^(x + e^, E, 9) - ^{x, E, 9)) 

With these definitions we arrive at the time evolution 

(17) E{t + l,x) = E{t,x)+V-3{x,E{t),9{t)) 

(18) 9{t + l,x) = g{9{t,x)) + h{x,9{t))). 

Note that flT7|l is a natural discrete space time version of ([5]). Let us discuss this 
iteration from a general perspective before making more specific assumptions of the 
perturbations. 

5. Fast Dynamics 

The iteration (fTS!) of the chaotic variables is autonomous. We shall assume the 
perturbation h is C^ with the following locality property 

(19) \deiy)h{x,9)\<ee~''\^~y\ 
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and Holder continuity property 

(20) \de(y)h{x,9)-deiy)h{x,9')\ < eJ^e-^d^-^l+l^-^Dl^l^) - e'(z)|. 

z 

These properties guarantee [8] that the 6'-dynaniics is space-time mixing. This means 
that the dynamics is defined in the A — )■ Z'' hmit and it has a unique Sinai- Ruelle-Bowen 
measure fi on the cyhnder sets of N^ which satisfies 

(21) E{F{e{t, x))G{e{0, x))) - E{F{e{t, y))EG{e{0, y)) < Ce-"(*+l^-^l) 

for Holder continuous functions F and G. Here E denotes expectation in fi. 

We conclude that sampling ^(0, ■) with n makes 6{t,x) random variables which are 
exponentially weakly correlated at distinct space time points. Therefore 6x{t) acts as 
a random environment ioi the slow variable dynamics f lT7|) . 

6. Quenched diffusion 

The previous discussion shows that we can view the current 3{x,E,6{t)) in the 
slow variable dynamics (TT7|) as a random field J(t, x, E) which is exponentially weakly 
correlated in space and time. We may thus rephrase the problem of deriving diffusion 
in deterministic dynamics as that of quenched diffusion in random dynamics. We want 
to show that the random dynamical system 

(22) E{t + l,x) = E{t,x) + V ■3{t,x,E{t)):=<^{t,x,E{t)) 

has a diffusive hydrodynamical limit almost surely with respect to the SRB measure /i. 
Let us inquire how this should come about and then list the assumptions we need for 
the actual proof. 

Consider first the annea/ec? problem, i.e. averaged equation (!22|) : 

E,{t + 1) - E,{t) = V ■ E[J(t, X, E{t))] := V ■ J{x, E{t)). 

where, by stationarity of /i, J' is time independent. Supposing that h and J have 
natural symmetries under lattice translations and rotations we infer that J^ vanishes 
at constant E and then locality assumptions of the type we assumed for h imply 

Jix,E) = J2<^^y^E)VEiy). 
y 
Hence the annealed dynamics is a discrete nonlinear diffusion 

Fit + 1) - Fit) = V ■ K{F{t))VF{t) 

provided the diffusion matrix n{F{t)) is positive. 
Let now 

(3{t, X, F{t)) = Jit, X, Fit)) - Jix, Fit)) 
be the fluctuating part. Then slow dynamics becomes 

Fit + 1) - Fit) = V ■ KiFit))VFit) + V ■ Pit, Fit)) 

E l3it,F) = 
i.e. a nonlinear diffusion with a random drift. In a physical model one would expect 
niFit)) to be positive although not necessarily uniformly in F. If furthermore /3 turned 
out to be a small perturbation quenched diffusion might be provable. In what follows 
we will make such assumptions and then indicate how to establish diffusion. 
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Before stating the assumptions let us make one more reduction. It is reasonable to 
assume E = is preserved by the dynamics. This then implies (3(t, 0) = 0. Let us 
study the linearization a.t E = 0: 

(23) E{t + 1) - E{t) = V ■ K{0)VE{t) + V ■ {D/3{0, t)E{t)) 
or, in other words 

(24) E,{t + l) = J2p^yit)Ey{t) 

y 
with 



J2p==yi^) = 1- 



Since E' > we have pxy > i.e. Pxy{t) are transition probabilities of a random walk. 
Pxy(t) is space and time dependent and random i.e. it defines a random walk in random 
environment. 

7. Random walk in nonlinear random environment 

Consider a random walk defined by the transition probability matrix Pxy{t) at time 
t. pit) = pit, oj) is a taken random defined on some probability space Vt. We suppose 
the law of p is invariant under translations in space and time. Define 

(25) ||E|| :=sup|E(a;)|(l + |a;|)'^+" 

X 

for some a > 0. Let, att = 0, ||£'||<oo. We say the walk defined by p is has a diffusive 
scaling limit if there exists C, k such that almost surely in uj 

(26) lim \\L'^E{Lh, L-) - Ct-^/^ El{- / yft\\ = 

where E*{x) = e~^ 1^^ . In other words 

L'^EiLh, Lx) ~ C't-'^'/'e""'/^"* 

as L — )• oo. 

We prove this for a non-linear perturbation of RWRE. Let us state the assumptions 
for the random dynamical system eq. fl22|) . We assume $ is C^ in \\E\\i < S and 
satisfies 

Positivity: $(E) > for E>Q. 

Conservation law: 

Y,Ht,x,E) = Y,Ex 



Weak nonlinearity 

.d'^^it.x.E). 

< e 
' (JEydE, ' 

Write the average map 



d^^{t,x,E)^ ^ „_|,_j,l_|,_,| 



E$(t, x,E) =Y^ T{x - y)Ey + o^E). 

V 

EUipticity: T generates a diffusive random walk on Z'^. 
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Write 

$(t, X, E) - E^{t, X, E) ■=V ■ b{t, X, E). 

Weak correlations. Assume 

(27) b{t,x,E)= Yl bA{t,x,E) 

A(ZZ''x[0,t] 

with 

and Ba, bs are independent if A fl i? = 0. 

Remark. A representation of the form fl27|) arises from the model we have discussed 
above with the proviso that b^, bs are independent only in the case the 9 dynamics 
is local, i.e. /i = in eq. (ITSl) . For the general h there is weak dependence that can be 
handled. 

Theorem 7.1. Under the above assumptions and 5, e small enough the random dy- 
namical system $i is diffusive, almost surely in uj. 

8. Renormalization group for random coupled maps 

The proof of Theorem 7. 1. [9] is based on a renormalization group method introduced 
in [To] and [11]. Let us introduce the scaling transformation Sl. 

(28) {SlE){x)=L'^E{Lx). 

where L > 1. Fix L and define, for each n G N, renormalized energies 

E^it) = SL,.E{L^^t). 
We can then rephrase the scaling limit fl26l) as 

lim U"^E{L^''t,L''x) = lim En{t,x). 

n— >oo 71— >oo 

En{t) inherits dynamics from E. We will call this the renormalized dynamics: 

E„(t + l) = $„(t,E„(t)). 
Explicitely we have 

$„(t) = SLn(<l>(L2"t + L^" - 1) o . . . o $(L2"t))^i-n. 
The dynamics changes with scale as 

$„+i = 7^$„ 

with 

7^$(t,■) = ^L$(tL0o...o$(tl)5^l 

with ti = L^t and ^^2 = L2(t + 1) - 1. 

7?. is the the Renormalization group flow in a space of random dynamical systems. 
We prove: almost surely the renormalized maps converge 

7^"/ ^ /* 

where the fixed point is nonrandom and linear. 

f*{E) = e^^E. 
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Moreover, the renormalized energies converge almost surely to the fixed point 

iiE„(t,-)-^i?:(-M)ii^o 

which is the diffusive scaling limit. 

These results may be summarized by saying that both the randomness and the 
nonlinearity are irrelevant in the RG sense. Let us finish by sketching the reasons for 
this. 

We start by considering the linear problem 

D^t,x,0)E = J2p-yit)Ey. 
y 
Then DR^ = p' with 

p'it),y = L\p{L\t + 1) - 1) . . .p{LH))L^Ly. 

Write 

P^y{t) = T{X -y)+Vy- C^y{t) 

with Ep = T and Ec = 0. Then, for p' = T' + Vc' we get 

(29) T\x -y)= L'^T^' {Lx -Ly)+r{x- y) 
where r is an expectation of a polynomial in c. For the noise we get 

(30) V.4^ = L^5^ 5^T*(Lx - u)VuCuv{t)T'''~''\v - Ly) + 7,,. 

t=l uv 

where 7 involves quadratic and higher order polynomials in c. 
Ignoring first r we get for the average fiow 

T„ = L"'^T^'"(L"-) 

f„(A;)=f^'"(A;/L"). 
Write f{k) = 1 - ck^ + o{k'^). Then as n -> 00: 

explaining the fixed point. 

Similarly, ignoring 7 the noise is driven by the linear map 

L2 

£c,,(0) = L'^-i Y. H ^*(^^ - u)cUt)T''"-'~\v - Ly). 

t=l uv 

The variance of Cc contracts: 

(31) E(£c)' ~ L"'E^. 

The intuitive reason behind this is the following. Take e.g. x = y = Q. For t of order 
L^, T^{Lx — m) ~ L-'^e-\^-'^/^_ Hence the u and the v sums are localized in an L cube 
at origin. Since c„„(t) has exponential decay in |-u — f | 

(32) £coo(0)~L-^-i5^5^c,„(t). 

t=l |m|<L 
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Since correlations of c decay exponentially in space and time ( 132|) is effectively a sum 
of L*^"*"^ independent random variables of variance L^^'*~^(Ec)^ thus leading to (!3T|) . 

Taking into account the corrections r and 7 in (!29l) and (!30l) we conclude that the 
variance contracts as 

E(c„)2 ~ e„ = L-"'^e. 

The iteration of the mean becomes 

(33) T„+i = L'^Tf(L-) + 0(e„). 

The fixed point is the same but the 0(e„) renormalizes the diffusion constant k at each 
iteration step (less and less as n — )■ 00). 

There is a problem however once we try to make this perturbative analysis rigorous. 
Deterministically the noise is relevant from (l32l) we see that ||£c||oo can be as big as 
0(L)||c||oo- This means that there are unlikely events in the environment where the 
random walk develops a drift. We write 

|c„(t,,E)|<L^"(-)-''^ 

Then Nn{x) can be (very) large, but with (very) small probability: 

Prob(iV„(x) > A^) < e"^^ 

with K large. 

Finally, to control the nonlinear contributions to $„ we show that the second deriv- 
ative D\(^ is irrelevant in all dimensions due to the scaling of E: 

7^<l>(t, X, E) = L'^{^{tL2) o ■ . ■ o $(ti))(Lx, L-'^E{-/L)). 

9. Towards Hamiltonian systems 

The coupled map lattices we have discussed are an alternative microscopic model 
with a local conservation law that under a macroscopic limit gives rise to diffusion. 
To be realistic they should however share some features with the Hamiltonian systems 
that are more familiar and physically relevant. From this point of view there is a lot 
missing from our analysis. 

The first problem to understand is to go beyond the perturbative analysis around 
E = (i.e. zero temperature). Then the equation i^^ picks also a driving term. 

The second unnatural assumption is the i?-independence of the 6 dynamics. In 
a realistic model rare configurations of E can slow down the 6 dynamics. Also the 
annealed system is probably not uniformly elliptic as we assumed and the random drift 
can create traps in the environment with long lifetimes. 

All these issues can and should be be studied with the renormalization group ap- 
proach sketched above. 
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